JOURNAL OF RESEARCH of the National Bureau of Standards— B. Mathematics and Mathematical Physics 

Vol. 65B, No. 3, July-September 1961 

Theory of an Accurate Intermediary Orbit for 

Satellite Astronomy 1 

John P. Vinti 

(May 16, 1961) 

This paper derives an accurate intermediary orbit of an artificial satellite of an oblate 
planet. The drag-free motion takes place under the action of a gravitational potential 
which fits the even zonal harmonics exactly through the second and approximately through 
the fourth, in the case of the earth. This potential leads to separability of the Hamilton- 
Jacobi equation. 

Two alternative sets of orbital elements are set forth. The first set is related directly 
to initial conditions, but requires numerical factoring of a certain quartic to evaluate some 
of the integrals. The second set, on the other hand, permits exact factoring of both quartics 
that appear, but is not related directly to initial conditions, so that its members can best 
be obtained by a least-square fit of the solution over many orbital revolutions. 

The final solution is given in terms of certain uniformizing variables, whose periodic 
terms are correct through the second order in the oblateness parameter and whose secular 
terms are exact, for the intermediary orbit. These exact solutions for the secular terms are 
expressed by means of certain rapidly converging series, with complete avoidance of elliptic 
integrals of the third kind. Sections 9 and 10 give a summary and a list of symbols. 

1. Introduction 

The author has introduced a gravitational potential [l] 2 for an axially symmetric planet 
which accounts rather accurately for its oblateness and yet still leads to separability of the 
problem of satellite motion. The resulting solution is expected to yield an "intermediary 
orbit" somewhat more accurate than those heretofore used, in that it accounts for all of the 
second zonal harmonic and for more than half of the fourth zona! harmonic. (Previous in- 
termediary orbits [2, 3] have accounted for only part of the second harmonic and have neglected 
the fourth harmonic.) The residual fourth harmonic, the odd harmonies, the tessera! harmonics, 
the lunar-solar forces, and aerodynamic and electromagnetic drag are then to be considered as 
producing perturbations of this intermediary orbit. 

Izsak [4] has already given an analytical solution for this intermediary orbit, with both 
periodic and secular terms correct through the second order in a certain oblateness parameter. 
His solution makes rather heavy demands on the reader's knowledge of linear fractional trans- 
formations and the theory of elliptic functions in the complex plane. The present paper 
avoids these complications, with elliptic integrals occurring only in the simple forms of the 
complete first and second kinds. Furthermore the resulting solution not only gives the periodic 
terms correctly to the second order, but gives the secular terms "exactly"; i.e., to arbitrarily 
high order. I wish to acknowledge very explicitly, however, that I am greatly indebted to 
Izsak for the introduction of one of the sets of orbital elements that I use. Knowledge of 
this set, which permits exact factoring of a certain refractory quartic, has influenced my treat- 
ment of the whole problem. 

2. Statement of Problem 

If p, 7], are the oblate spheroidal coordinates introduced in [1] and if r, 6, $ and X, Y, Z 
are the corresponding spherical and rectangular coordinates, then 

X+iY=r cos 6 exp i<j> = [(p 2 +c 2 )(l-ri 2 )]i exp icf>, (2.1) 

Z=r sin d = PV) (-l^r/^1). (2.2) 



1 This work was supported by the U.S. Air Force, through the Office of Scientific Research of the Air Research and Development Command. 

2 Figures in brackets indicate the literature references at the end of the paper. 
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Here r is the geocentric distance to the satellite and 6 and <$> are respectively its geocentric 
declination and right ascension. For sufficiently large r, p^r and v\ — sin 6. 

With the origin taken at the planet's center of mass, the intermediary orbit is then the path 
of a particle in the approximate potential field 

V.—mpCp'+cV)" 1 , (2.3) 

where /* is the product of the gravitational constant and the planetary mass and where 

c 2 = rlJ 2 . (2.4) 

Here r e is the equatorial radius of the planet and J 2 is the coefficient of the second zonal harmonic 
in the expansion of the planet's true potential in spherical harmonics. For the earth J 2 — (1.08) 
10~ 3 ; to three significant figures. 

According to [1], the coordinates p, rj, </> satisfy the following equations, involving quad- 
ratures : 

t+fa=± P p*F-\p)dp±<? P v 2 G~Kv)dv, (2.5) 

J p 1 Jo 

ft= Ta 2 P F'Kf>)dp±a 2 f" 6-i(v)dv, (2.6) 

J p x Jo 

4>-/3 3 =Tc 2 a 3 P (p i +c*)- l F-i(p)d P ±a 3 P (1-r, 2 )" 1 ^-*^)^- (2.7) 

J P { Jo 

Here F(p) and G(rj) are the quartic polynomials 

F(p) = c*al+(p*+c*)(-a 2 2 +2w+2a l p*), (2.8) 

G( V ) = -<4+(l-7 i *)(<4+2a 1 c*r J *). (2.9) 

The as and fts are the Jacobi constants, with the energy a^O for satellite motion and with the 
polar component of angular momentum a 3 ^0 according as the orbit is direct or retrograde. 
To orient oneself, note that in the limiting case c— >0 of Keplerian motion the separation constant 
a 2 reduces to the total angular momentum, —ft to the time of passage through perigee, ft to 
the argument co of perigee, and ft to the right ascension 12 of the ascending node. The Jacobi 
constants may be determined, at least in principle, from the initial conditions; we have more 
to say about this point later. The constant pi is the next-to-the-largest real zero of F(p) and 
thus is that zero of F(p) which is closest to the smaller zero of 

f( P )^-a 2 2 +2^p+2a 1 p\ (2.10) 

To solve (2.5) through (2.9) for p, 77, and 4> as functions of t, we must first solve (2.5) and 
(2.6) for p and rj and then substitute the results p(t) and rj(t) into (2.7) to determine <£(£). To do 
so we must first evaluate the above six integrals, which we shall obtain in terms of certain uni- 
formizing variables. In turn, evaluating these integrals presupposes knowing how to factor 
the quartics F(p) and G(rj). 

3. Factoring the Quartics: Orbital Elements a , e , io, ft, ft, ft 

In the case of elliptic motion (c=0) the perigee and apogee radii r x and r 2 would be the two 
zeros oij(p), viz, 

n=a (l— e ), (3.1) 

r 2 =a (l+eo), (3.2) 
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where 

ao^-i/^f 1 (3.3) 

6^1+2^ 2 V 2 . (3.4) 

In our present problem, with c^O, we may still define constants a and e by (3.3) and 
(3.4), as well as another constant 

^eecos" 1 (ajj/afc). (3.5) 

The constants a , e , i , ft, ft, and ft constitute one ijossible set of orbital constants. We may also 
introduce the corresponding semi-latus rectum p , denned by 

Po = a (l-4) 9 (3.6) 

so that 

<4=m- (3-7) 

A determination of a ]} a 2 , and a 3 would then furnish a , e , and i . If the subscript i denotes 
initial values, then (per unit mass) 

ai=ltff-w«(p?+A?J)- 1 , (3.8) 

a z =r\ cos 2 Bd^Xtti—YiXt, (3.9) 

where u means speed and a superscript dot denotes the time derivative. Also, by equations 

(50), (59.1), (13.2), (10.2), and (55) of [1], 

aH(i-^)- 1 [(p!+^) 2 ^+^-2a lC 2 ^(i-^)]. (3-10) 

Thus a knowledge of the initial coordinates and their initial derivatives would determine the 
as and thus the constants a , e , and i Q . 

A knowledge of their numerical values would then permit a numerical solution of the 
quartic equation F{p)~0 and thus furnish the numerical values of p u p 2 , A, and B necessary to 
factor F{p) into the form 

F(p)^-2a l (p- Pl )( P2 - P )(p*+Ap + B), (3.11) 

where pi and p 2 are the zeros of F(p) closest to the values r x and r 2 . Then, in the intermediary 
orbit, p is restricted to the interval pi ^ p ^ p 2 between two spheroids. 

By equating the coefficients of corresponding powers of p in (3.11) and (2.8) we find 

P 3 : p l +p 2 -A=- f jLai l =2a , (3.12) 

p 2 : B+p 1 p 2 -(p l +p 2 )A=c 2 -i<4c*r 1 =C*+a po, (3.13) 

p: (p 1 +p 2 )7i-p lP2 /l=- M cV 1 =2a c 2 , (3.14) 

p°: pip 2 B= — %c 2 (al—al)a{ 1 =a poC 2 sm 2 i , (3.15) 

with use of (3.3) through (3.7). 

By beginning with the zero-order solution A= B=0, Pi J rP2=r 1 J r r 2 =2a , p\P2=r l r 2 = a {) pQ J 
one can solve this set of equations for the four unknowns A, B, P1 + P2, and pip 2 , by a 
method of successive approximations. If 

h^c 2 /pl^(r e /p ) 2 J 2 , (3.16) 

x=(l-et)$, (3.17) 

y=a 3 /a 2 =cos i , (3.18) 
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the second-order solution, through kl, is 

A=-2k p y 2 [l+k (2x 2 -Sx 2 y 2 -4:+8y 2 ) + . . .], (3.19) 

B=k p 2 (l-y 2 )[l+h(4y 2 -x 2 y 2 ) + . . .], (3.20) 

Pl +p 2 =2p Q x~ 2 [l-k x 2 y 2 -k 2 x 2 y 2 (2x 2 -3x 2 y 2 -4,+Sy 2 ) + - - .], (3.21) 

PiP2=^~ 2 [l+% 2 (^ 2 -4)-^(12x 2 -x 4 -20^V-16+32^ 2 +xy)+ . . .]. (3.22) 

The constants 

as*(pi+Pa), (3.23) 

e=^^, (3.24) 
P2 + P1 

^ = a(l-6 2 ) (3.25) 

will occur throughout the evaluation of the p-integrals. In terms of a , e , and i , their values, 
to the second order in k , then satisfy 

^ =l+2&o2/ 2 (x 2 -2)+% 2 (3^ . . ., (3.26) 

p /p=l+2k y 2 (2-x 2 )+kly 2 (-3x*+6xy+mx 2 -40xY-W^ . . ., (3.27) 

1-e 2 



2\1 



--l+koy 2 (3x 2 -4)+k 2 oy 2 (5x*-2xy-20x 2 +28x 2 y 2 +lQ-32y 2 )-{- . . ., (3.28) 



(^^Y = 144^ 2 (3x 2 -4)+^ 
\1 6q/ £ 8 



(3.29) 



(l-^ 2 )^- 1 =a - 1 (l-^)-^[l+i/:oy 2 (4-x 2 )+iA:y(-4x 4 +7xV+48x 2 

-104xV-64+176y 2 )+ . . .]. (3.30) 

In those places where e occurs alone in the theory, i.e., not in the combination 1 — e 2 , it may 
of course be found by use of 

e = [l-(l-e 2 )]K ( 3.31) 

By (3.28) this results iu 

e=[e 2 +k x 2 y 2 (4:— 3x 2 )+ . . .]*, (3.32) 

so that when e is comparable to k it is not feasible to expand e in a power series in k 'j indeed 
if e =Q, we should need a power series in kl. 

Direct use of this second-order solution in factoring F(p) will lead to p-integrals that have 
secular terms correct only to 0(kl). Since we are aiming at arbitrarily high accuracy for the 
secular terms, we include it here for other purposes. The most important of these purposes 
is to furnish information about the orders of various quantities in k ; e.g., A and B are both of 
order k . Such information will be necessary in carrying out the solution of (2.5) through (2.7) 
for the periodic terms. The second purpose is to furnish a convenient starting point for any 
investigators who may choose to use a , e , and i , along with the /3's, as orbital elements, and 
who will therefore have to solve the equation F(p) = Q numerically. A third purpose is for use 
in calculating the mean motions to first order, for comparison with other theories. 

The quartic G(r}), which is quadratic in rj 2 , may be factored either as 

G^^^a^ivl-^irjt-rj 2 ) (3.33) 
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or as 

G( V )^(al-adv 4 (v- 2 -Vo 2 )(v~ 2 -V2 2 ). (3-34) 

If we use the latter form, we find on comparing (3.34) with (2.9) that rjQ 2 and ^ 2 ~ 2 are the roots 
of the following quadratic equation in r\~ 2 : 

(a 2 -a 3 2 )?T 4 + (2a 1 c 2 -a 2 2 )r 1 - 2 -2a 1 c 2 =0. (3.35) 

Thus 

= J(«2-2aiC 2 )(^-ai)- 1 (l+^), (3.36) 

4(^-2a 1 c 2 )(^-al)- 1 (l-^), (3.37) 

where 

W^l + 8a 1 c 2 (a 2 -ai)(a 2 -2a 1 c 2 )- 2 . (3.38) 

From (3.36) and (3.38) it follows that for ai <0 

Vl^^ 2 (a 2 -ai)^\ (3.39) 

and we shall see below that for satellite motion i^^l. Since yf cannot exceed 1, it follows that 
in the actual motion rj always lies in the interval — 770^ v^ Vo between two hyperboloids. 
We readily find that 

77o=(sin io) [l-hhxy+iWxy(7y*-±)+ . . .], (3.40) 

(1 — ^c 2 ))-^=|sec^|[l — i^ x 2 (l— 7/ 2 ) + |^xXl— y 2 )(57/ 2 — 1)+ . . .], (3.41) 

v ^=koX*(l-k xy +...), (3.42) 

(V^) 2 =*^ 2 (sin 2 io)(l-2k xy+ . . .). (3.43) 

Note that r)2 2 ^k ^10-\ so that 77;^ 1000. 

4. Factoring the Quartics: Orbital Elements a, e, I, ft, ft, ft 

If we equate the coefficients of powers of rf in (3.33) with those of corresponding powers 
in (2.9), we find 

vl+vl=l~^ 2 =l+a Po/c*, (4.1) 

2 2 

If in (3.12) through (3.15) we use (3.23) through (3.25), we find 

2a-A=2a , (4.3) 

B+ap— 2Aa=c 2J ra p Q , (4.4) 

2aB-Aap=2a c c 2 , (4.5) 

Bap=a p c 2 sm 2 i . (4.6) 

Suppose we now regard a, e, and 

77o = sin / (4.7) 

as known. (When we adopt a, e, I, and the /3's as orbital elements we are certainly assuming 
so; we discuss later how they may be determined.) Then in treating the p-integrals we have 
five unknowns, viz, a , e , i , A, and B, and in treating the 17-integrals one additional unknown, 
viz, 7] 2 . Altogether then, we have six unknowns and we have six equations with which to 
determine them, (4.1) through (4.6). With these orbital elements, viz, a, e, and rj , however, it 
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turns out that the system can be solved exactly and with considerable ease, for the required 
unknowns. This property of a, e, and % was first pointed out by Izsak [4]. 

To carry out the solution, first eliminate -q 2 from (4.1) and (4.2). The result is 

?I^Lsin* i = v l+ cW - vl) - (4.8) 

a 2 a p 

On inserting (4.8) into (4.6), we find 

Bap=cS\ [aoPo+c 2 (l-ri)). (4.9) 

If we now use (4.3) to eliminate a from (4.5) and (4.9) to eliminate a p from (4.4), we find a 
pair of simultaneous linear equations for A and B: 

(ap-c 2 )A-2aB= -2ac 2 , (4.10) 

2c 2 r ] laA+(ap-c 2 rjl)B=c 2 v 2 (ap-c 2 r 1 2 ). (4.11) 



Their solution is 



Then, from (4.3) and (4.12) 



A= 2ac 2 (]-r } 2 ) )(ap-c 2 ri 2 ) 12) 
(ap—c 2 ) (ap—c 2 rjl)+4a 2 c 2 r]f 

„ . 2 (ap-c 2 )(ap-c 2 v 2 )+4a 2 c 2 , 

Vo (ap-c 2 ) (ap-c 2 rj 2 )+4a 2 c 2 v 2 ( *' l6) 



ix _a A c 2 (l-rjl)(ap-c 2 r]l) . 

~ „ 1 Or, L ' /™ „2\f„^ „2„2\\A„2„2»2 K^- 1 -*) 



2aai a 2a (ap—c 2 )(ap—c 2 rj 2 )+4a 2 c 2 7}l 

and from (4.9) and (4.13) 

Equation (4.8) and the relation tj = sin I give 

( c 2 vl\* 
0Lz=a 2 \l — -^pj cos /, (4.15a) 

so that (4.15) and (4.15a) determine a 3 . Finally, to obtain rj 2 , combine (4.2) and (4.6) to 
obtain rj2 2 =c A 7jl (Bap)' 1 and then use (4.13). The result is 

_ 2= j^ (ap-c 2 ) (ap-c 2 r) 2 o)+4:a 2 c 2 >n 2 

2 ap (ap—c 2 )(ap—c 2 r)D J r 4a 2 c 2 ^ 

This completes the solution for the required unknowns when the orbital elements are a, e, 
and /. In terms of these orbital elements we can now factor the two quartics F(p) and G(rj) 
exactly and thus evaluate all the integrals. 

With use of the old oblateness parameter k = c 2 /pl and a new one, suitable for use with this 
second set of orbital elements, viz, 

Jc=c 2 /p 2 ^tr e /p) 2 J 2j (4.17) 

we can readily show that, to the first order, the equations of sections 3 and 4 give similar results. 
Thus we readily obtain 

A** — 2k p cos 2 i « — 2kp cos 2 I, (4.1 8) 

B~hpl sin 2 i Q ~kp 2 sin 2 /, (4.19) 

-«l+Jfco(l-4>) cos 2 ^l+£(l-6 2 ) cos 2 /, (4.20) 

a 
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7 -^«14-2& (l+« 2 ) cos 2 i ~l+2k(l + e*) cos 2 I, (4.21) 



l-el 
1-e 2 



'-l+k (l+3el) cos 2 i «l+fc(l+3e 2 ) cos 2 /, (4.22) 



^5?«l+k,(l-«J) cos 2 i «l+^(l-e 2 ) cos 2 7, (4.23) 

Sill 1 jLl & 

v7~h(l-e&+ . . . ~£(l-6 2 )+ . . , (4.24) 

(Vo/ri2) 2 ~ko(l-el) sin 2 i ~&(l-e 2 ) sin 2 /+ (4.25) 

With either set of orbital elements, the final solution will be given in terms of a, e, I, ft, 
ft, and ft, and certain angles E, v, and ^, analogous respectively to the eccentric anomaly, the 
true anomaly, and the argument of latitude in elliptic motion. Once one knows a, e, and I, 
one can then determine the fts by observing, at various times, whatever quantities will best 
serve to determine E, v, and \f/. 

If one is using a , e , and i as elements, one has to determine them from initial values or 
from some procedure equivalent to determining initial values. One then lias to factor the 
quartic F(p) numerically to del en nine a, e, and 770 — sin /. 

if one is using a, e, and I as elements, one has to determine them by following the orbit 
for many revolutions and then applying some sort of iterated least-square process. In this 
case one can then find the Jacobi constants ot\, a 2 , and a 3 and thus a , e , and i by means of the 
equations of this section. 

In any event, the determination of orbital elements, by comparison of theory and ob- 
servation, is ordinarily considered a completely separate problem in celestial mechanics from 
the calculation of the motion for given orbital elements. We have included the above remarks 
only to aid in the possible application of the present theoretical solution. Indeed the problem 
is further complicated by the small perturbations that occur in practice. 

In the rest of the paper we simply assume a, e, 77 = sin 7, ft, ft, and ft to be known and t hen 
complete the solution for the intermediary orbit. There will be one restriction, however. 
The method of evaluating the p-integrals will be found to depend on the orbital inclination 
?'o or /. For values of I less than 1°54 / or greater than 178°6 / , a different approach would be 
needed. We shall therefore restrict considerations in the present paper to orbits that have 
inclinations between these two values, thereby ruling out equatorial or almost equatorial 
orbits. 

5. The p-Integrals 

In (2.5) through (2.7) the p-integrals are 

Ri = ± ( P P 2 F~Kp)dp, (5.1) 

R 2 = ± ( P F~Kp)d P , (5.2) 

Rs = ± ( P (p 2 +c')- 1 F-Hp)dp, (5.3) 

where F(p) is given by (3.1 1), A and B by (4.12) and (4.13), p by (3.25), and p, and p 2 by 

Pi=a(l — e) f p 2 =a(l+e). (5.4) 

Equations (5.4) follow from (3.23) and (3.24). Then 

F-Hp) = (-2a l )-h(p-pO(P2-p)}-h- 1 (l+Ap-'+Bp-')-K (5.5) 
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If, for convenience, we put 

Am-2b u B = b 2 23 X^bi/h (5.6) 

and 

h = b 2 / P , (5.7) 

then 

(l+A P ~ 1 +Bp- 2 )-^(l-2\h+h 2 )-^ (5.8) 

=Jbh n P n (\), (5.9) 

if |A|<1 and |X|<1. 

The expansion (5.9) will be a convenient tool for evaluating all the p-integrals, whenever 
it can be used. To see when, use (5.6), (5.7), (4.18), and (4.19). Then to 0(k) 

bi=kp cos 2 /, b 2 =k*p sin /, 
so that b 2 /p<^.l and 

X=&* cos 2 / esc I^>0. 
Thus to 0(k), |X|<ll whenever 

k$ cos 2 / esc 7<0, 
or whenever 

tan 2 /+tan 4 />A: 
or 

2 tan 2 />(l+4£)*-l. 
To 0(k) } this becomes 

tan 2 I>k. 

But k=(r e /p) 2 J 2 and e/ 2 =(1.08)10~ 3 for the earth, so that \<1 whenever 

|tan/|>0.033r e /p. 

For close orbits r e /p~l, so that the necessary condition becomes 

I C <I<1S0°-I C , J c «l°54'. 

The expansion will thus work for all satellite orbits that are inclined more than 1°54 / to the 
equator. For those orbits that lie closer to the equator one must use some other method to 
evaluate the p-integrals. On inserting (5.7) and (5.9) into (5.5), we find 

F(p)-*=(-2oi)-* 2 &SP.(X)p" 1 -[(p-pi)(p2-p)]"*. (5.10) 

With use of (5.1), (5.10), and (5.6), we find 

(-2fl*)*Si= f (p+W[(p-Pi)(P2-p)]"*(±dp)+S &2 n P»(X) P P 1 - n [(p-pO(P2- P )]-K±dp). 

J Pl n = 2 J p l 

(5.11) 

Each of the separate integrals in (5.11) is a multiple-valued function of p. It is appropriate 
to change variables in each to a uniformizing variable E or v, analogous respectively to the 
eccentric and true anomalies in elliptic motion. We define E and v by requiring them to 
satisfy 

p=a(l — e cos E) = (l J re cos v)~ l p (5.12) 

and always to increase with time. Then, from (5.12), (5.4), and (3.23) through (3.25), we find 

[(p-p 1 )(p 2 -p)]-i(±dp)=dE=(l-e 2 )*(l+e cos v)~ l dv (5.13) 
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5.1. The Integral Ri 

On introducing E into the first integral in (5.11) and v into each of the integrals in the 
series, we find 

(-2a 1 ^R 1 =b 1 E+a(E-esmE) + (l-e 2 ^pJ^ (hlp) n P n W (\l + ecosv) n -*dv. (5.14) 

n=2 Jo 

To investigate the convergence of the series we write 

Si = J^ (b 2 lp) n P n {\) (" (l + e cosv) n -*dv (5.15) 

71 = 2 J0 

= (hlp) 2 Z) (hlp) m P m+ 2(\){ ' (1 + e cos v) m dv. (1.16) 

m=0 J 



Then 



\Si]£(h/p)*j: (b 2 /pr\P n+2 (\)\(l+e) m v (5.17) 

w=0 



and since |P n (X)|^ 1 for all |X|^ 1, we have 

\Si\£(kJp)H>±[!&££T (5.18) 

ro-o L P J 

< (^) 2g . (5 19) 

By (4.19) and (5.6), however, b 2 /p~k% sin /and of course 1 + e ^2, so that 

& 2 (l+e)/2>g2£* sin 7g2R (5.20) 

Since & — 10 -3 , we have 62(1 +^)//> ^0.063. Thus the series Si converges, and converges more 
rapidly than a geometric series of common ratio ~ 1/16. Actually, since we have shown that the 
series Si converges absolutely, we can regroup it into the sum of a series S ie containing only the 
even values of n and a series S\ containing only the odd values of n. It is then a simple matter 
to show that Si = S u -\-Si converges more rapidly than the geometric series [l-\-b 2 (l+e)fp]. 
2»!o ^2(1 + e )/p] 2n - That is, we can actually expect the convergence of S { to be as rapid 
as that of a geometric series of common ratio [b 2 (l J re)/p] ^4&^ 1/250. 

To decompose the series Si into a part proportional to v and a periodic part, note first that if 

Mv)=f\l+e cosv) m dv, (5.21) 

then f n (v) — (2Tr)~\f m {2w)v is an odd function of v, of period 2w. Hut f m (2ir)=2f m (w), so that 
we obtain 

fv fir 

f m (v) = \ (l+e cos v) m dv=w- 1 v (l+e cos v) m dv+^c mj sin jv. (5.22) 

Jo Jo j= 1 

To obtain, to any order in k, the parts of the p-integrals proportional to v, we shall need to con- 
sider all integral values of m in 7r _1 / m (7r). To obtain the periodic parts correct to order k 2 y we 
shall need values of m only up to 4. 

To obtain a convenient expression for the ^-term o(f m (v), note that [5] 



/; 



(z+^z 2 -l cos v) m dv=irP m (z), (5.23) 



for all values of z, including real values greater than unity. Here P m (z) is the Legend re 
polynomial 

p ^lrt^- l)m > (5 - 24) 

59668G— 61 2 177 



of the same polynomial form in z that holds when [2|<0, when it can be defined by the usual 
generating function. If we put 

g=(l-e")-* (5.25) 

in (5.23), we find 



i 



where 



(l + e cos v) m dv=T(l-e 2 ) m/2 P m [(l-e 2 )-*] 
R m (x)^x m P m (lM (O^x^l). 



(5.26) 
(5.27) 
(5.28) 



Thus R m (x) is a polynomial of degree [ra/2] in x 2 . The first few of these polynomials are given 
in table 1. 

Table 1 



R m (x) 



1 

H(3-j2) 

3^(5-3*2) 

(H) (35-30^4-32:4) 

( 1 / i)(63-70j24-i5z4) 

(He) (231 -315.t2+ 105^4—5^) 



From (5.22) and (5.27) we then obtain 



r*v co 

(l + e cos v) m dv=R m (-Jl — e 2 )v+^, c mj • sin jv. 

Jo ;=1 



(5.29) 



Through m=4, the coefficients c mj are easy to find, simply by binomial expansion of the integrand 
and conversion to a trigonometric polynomial. 
The results are given in the following table: 



Table 2. ( 


Coefficients c m i in {5.29) 


V 

m \ 
\ 


1 


2 


3 


4 



1 
2 
3 
4 


1 

2e 
3e-(-3e3/4 
4e+3e3 


e2/4 

3*2/4 

3e 2 /2+6V4 


<?3/12 

eV3 


eV32 



In (5.14), b n 2 = 0(k n/2 ) and P n (\) = 0(\°) = 0(k°) if n is even or P w (X) = 0(X) = 0(W) if n is 
odd. On inserting (5.29) into (5.14), using table 2, and keeping periodic terms through 0(k 2 ) 
only, we find 

(-2a 1 )^R 1 =b l E+a(E-e sin E)+A 1 v+^A lj -smjv, 



where 



and 



;'=i 



A 1 =(l-e^p± l (b 2 /p) n P n (\)B n ^l--e i ) 



^n=^^(-2b l b 2 p+bi)e } 



^ 12 ~ 32p 3 6s6 ' 



(5.30) 

(5.31) 

(5.32) 
(5.33) 



The above proof of convergence of the series for R x also shows the rapid convergence of the 
series (5.31) for the coefficient A x of the ?;-term. 
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5.2. The Integral R> 

On inserting (5.10) into (5.2) and using (5.12) and (5.13), we find 

(-2a 1 )7i > 2 =(l-6 2 )V 1 iZ(b 2 tp) n P n (\) C (1 + e cosv) n dv. (5.34) 

n=0 Jo 

As above, one can show at once that this series converges more rapidly than the geometric series 
fir=oWl+^)/ri n , with common ratio b 2 (l J r e)/p 5^2^-0.066. The same proof then applies 
to the coefficient A 2 in 

4 

(-2a l yR 2 =A 2 v+^A 2J smjv. (5.35) 

7 = 1 

Using the same methods as for R ly we then find 

A 2 =(l-e*yp- l I2(b2/p) n P n (\)R n (^-e*) (5.36) 

B = 

and, through periodic terms of 0(k 2 ): 

A 21 =(l-e 2 ) i p- 1 e[b,p- l + (3bi-b'f)ir^-^bl(\ + ey4)2r i +l-bi(4:+3e i )p- i l (5.37) 

^2 2 =(l-e 2 )V 1 [^(36T-^)/r 2 -§^6 ! 61/r'+|^(6e 2 +e 4 )/r 4 ]- (5.38) 

^2 3 =(l-e 2 )^-' ~ i-bMv-'+Kir'), (5.39) 

O 

A 2i =(l-e*y*p-*~bie\ (5.40) 

5.3. The Integral i? 3 
If in (5.3) we now use the binomial expansion 

(p'+c 2 )- 1 ^- 2 f] (-l) j c 2j p- 2j (5.41) 

;=o 

and insert the expressions (5.10), (5.12), and (5.13), we find 

(-2oO*i?3=(l-e 2 )^ 3 r;S(-l^^^ (5.42) 

Jo ;=0 n=o 

where the integrand is the product of two series, each of which converges absolutely for any 
value of v. Then [6] it is equal to the series formed by summing the products of the individual 
terms, taken in any order, and this resulting series is itself absolutely convergent, for any value 
of v. It is therefore uniformly convergent, by the Weierstress M-test [7], so that it may be 
integrated term by term. 

Let us now rewrite (5.42) in the form 

(-2a 1 )^ 3 =(l-6 2 )^- 3 P JC D n (l + e cos v) m + 2 dv. (5.43) 

Jo ra=-0 

Here 

D m =J}djd n ', (5.44) 

summed over all those non-negative integral values of j and n' for which 

2j+n'=m, (5.45) 
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and where 

dj={-iy{cipy, (5.46) 

a«'=(V2>) n 'iVP0- (5-47) 

Then, because of the uniform convergence, 

(-2ai)*#3=(l-e 2 )*p- 3 Zl Dm f" (1 + 6 cos ?;) w + 2 <fo. (5.48) 

TO=0 Jo 

It is interesting to investigate here the rapidity of convergence of the series 

S= S Dm f ' (1 + « cos ») ra+2 *>. (5.49) 

71=0 Jo 

If m is even, we have m = 2i and n' = 2n, so that 2j+2^ = 2i and j=i—n. Then 

Z>*=A«=S*-.Ai.=:E^ (5.50) 

ti=0 71=0 

so that 

I Ai|£ (<#) 2 ' S (& 2 /c) 2n . (5.51) 

71 = 

But c 2 =kp 2 and bl^kp 2 sin 2 7, by (4.19), so that 

(6 2 /c) 2 «sin 2 Jgl. 
Then 

|A«|^(»+l)fc'- (5.52) 

If m is odd, we have m = 2i+l and n' = 2n-\-l, so that 

D 2i+1 =i d i . n 8 2n+1 =± l (-iy- n (c/p) 2i - 2n (b 2 / P y^P 2n+l (X). (5.53) 

72=0 72 = 

Then 

|A«+i|£ (c/p) 2t (h/p) S (,b 2 /cy^k'ki (i+l) sin/. (5.54) 

72 = 

Then, breaking up S into an even series £ e and an odd series S , we find 

OO 00 

l&I^S (i + l)^(l + e) 2i+2 ^(l + <0 2 z;;£ (i+l)[jfc(l + e)*]'. (5.55) 

i=0 t=0 



Using S;x i =(l-j;)- 1 and Sru'=a;(l-a!)- 2 , we find So"(«'+l)a!-(l-^" ! , so that 

|o|<_(l±iO^_ 

1 e| = [l-£(l+e 2 )] 2 " 
Similarly 



^(l + e) 3 ysin i 
= [l-&(l+e) 2 ] 2 
Thus, 



(5.56) 



\S,\^k*(l+e) a va.n / S (i+l)[fc(l + c)*]' (5.57) 

1 = 



(5.58) 



I <?l < (l+e)Ml+^ (1 + e) sin/] 

The series (5.49) for R z thus converges more rapidly than the series expansion of the function 
given in (5.59). 
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Through 0(k 2 ) the values of D m are given by 

Table 3 



m 


Order 


z>™ 



1 
2 
3 
4 


ft 
k 

fc2 


c?o5o= 1 

dodi = 5i=bi!p^k cos 2 7 
do5 2 +rfi«o=(62/p) 2 P2(6i/62)-(c/p)2 
do53+di«i = (&2/p) 2 P3(V&2)-(c/p)2(6i/P) 

do64+rfi62+rf25o=(6 2 /p)^4(6i/62)-(c/p)2(6 2 /p) 2 P2(6,/&2) + (c/p)< 



As with i?x and i? 2 , we find 



where 



(-2a,)^ 3 =^ , +S 4u sin>, 



^ 3 =(l- e 2)^- 3 T. D m R m+2 (^l-e*). 



(5.60) 

(5.61) 



The remarks about the convergence of the series (5.49) appty also to (5.61). 

Since R 3 is multiplied by c 2 =&y; 2 in (2.7), we need periodic terms only througli 0(k), in 
order to have periodic terms in the final solution correct through 0(k 2 ). By (5.43) and tables 
2 and 3, their coefficients are 



A n =(l-e*)lp- z e [2+6,/r 1 (z+^-p-* (|+c 2 ) (4+3e 2 )], 
A 32 = (l-e>y>p-> g+| h e*-p^+p (f +* 2 )} 

6. The Tj-Integrals 
In (2.5) througli (2.7) the ^-integrals are 

N 2 =± £ G-Hv)dn, 



(5.62) 
(5.63) 
(5.64) 
(5.65) 



N 3 =± £ (l-v*)- l G(t,)-*d v , 



(6.1) 
(6.2) 
(6.3) 



where G{r,) is given by (3.34), % by (4.16), a\ by (4.14) and (4.15), and a\-a\ finally by (4.8). 
In evaluating Ni and N 2 it is convenient to put 



tj = ))o sin rp, 



(6.4) 



where ^is to be an angle that always increases with time. In the limiting case c=0 we should 
have sin ^=sin0/sin / and ^ would thus reduce to the argument of latitude, i.e., to the angle 
between the line of nodes and the radius vector to the satellite. 
Then 

±G-Hv)d V =(4-al)-ho(l-q 2 sm'4')-W, (6.5) 
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where 

g=Wla=0(**)<l, (6.6) 

7j 2 being given by (4.16). 

6.1. The Integrals N x and N 2 
Insertion of (6.4) and (6.5) into (6.1) gives 

Ni=(c$-ci)-*fil fV-2 2 sin 2 $"* sin 2 ftfy. (6.7) 

Because of the identity 

(1 — q 2 sin 2 yf)~^q 2 sin 2 ^=(1 — q 2 sin 2 \p)~^— (1 — # 2 sin 2 $*, (6.8) 

M=(«i-o|)-^~ 2 [" ( V~2 2 sin 2 £T*#- fV-2 2 sin 2 *)*#| (6.9) 

= (<*!-<*§) -Vo<T 2 [F^,q)-E(M)}, (6.10) 

where F(\//,q) and E(\f/,q) are respectively the elliptic integrals of the first and second kinds, with 
modulus <£. 

Insertion of (6.4) and (6.5) into (6.2) gives 

N 2 =(al-4)- h Voj\l-q 2 sin 2 ^)~*# (6.11) 

= (c 2 -a 3 2 )-^ F(^g)- (6.12) 

Thus iVi and iV 2 can both be expressed exactly in terms of elliptic integrals of the first and 
second kinds. Our ultimate purpose, however, is to express each uniformising variable as the 
sum of an exact secular term and periodic terms correct through 0(k 2 ). For this purpose it is 
convenient to express each elliptic integral as a linear combination of \p and a Fourier series 
SrB» sin 2ra*. 

To do so, consider 

is I (l-q 2 s'm 2 x)~>dx. (6.13) 



F(M) ^£ 



Some simple transformations show that 

FU,+T,q) = FU,,a)+2K(q), (6.14) 

where 

K(q)= J ' \l-q 2 sin 2 x)~Hx (6.15) 

is the complete elliptic integral of the first kind. It follows that the function F(\(/,q) 
— (2lir)K(q)\p is periodic in \f/ with period x. 

Furthermore it is an odd function of ^, so that it can be expanded in a Fourier series containing 
only the sines of even multiples of \f/. Thus, 

F{^q) = {2lTr)K{q)t+itF qm sin 2m* (6.16) 

To calculate the Fourier coefficients F qm , differentiate (6.16) with respect to ^ and use (6.13). 
Then 

(l-2 2 sin 2 0-»=(2/x)#(2)+2Jb»»/ ? Ucos2m^. (6.17) 

i 
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The Fourier coefficients F gm are then given by 

F Qm =(2/wm) (1-g 2 sin 2 z)-* cos 2mzcte. (6.18) 

Since cos 2ra# is a polynomial in cos 2 x, each F Qm can ultimately be expressed as a linear com- 
bination of K(q) and the complete elliptic integral 



E 



?(g)= (l- 2 2 sin 2 x)^x. (6.19) 

Such a procedure, however, would not readily reveal the order of each coefficient in q 2 , which 
itself is of order k. (Actually we shall show that F Qm is of order k m , so that we shall need only 
F a i and F q2 .) Instead we expand the radical, obtaining 



Then 



Also 



and 



(l-g 2 sm 2 z) *=l+g- — 2 2n (n\) 2 ( ^ 

2 °° (9flVfl 2n fx/2 

i^m= — S on/' im sin2w x cos 2mz(/jc. (6.21) 

sin2W ^ = (§^ +( ~ 1)n21 " 2n § ( ~ 1} ' (2^— i)!i! C0S ^~ 2 ^^ ^=^ < 6 ' 22) 

cos [(2n—2j)x\ cos 2raz=| cos [(2n+2m— 2j)x]+% cos [(2^— 2m— 2j>]. (6.23) 



The integral of this product from to t/2 fails to vanish only if j=n-\-m or if j=n—m. The 
value j=n-\-m is absent from (6.22), so that only the term j=n— m in (6.22) contributes to 
the integral in (6.21). Thus 

( ' 2 sin 2w x cos 2mxdx=(-l) m 2- 2n ^ — rr £, (n>m). (6.24) 

Jo (n+?n)\(n—m)\ 2 v " 

On inserting (6.24) into (6.21), we find 

v 7 ^T ft (^+m)!(n— m)!2 4w (w!) 2 

Inspection of (6.25) now shows that F qm =0(q 2m ) = 0(k m ). The first two values are, to order 
q l or order k 2 : 

^ = -|(l+j5 2 )+.--, (6.26) 

F *=i£ 6 +- • • • ( 6 - 27 > 

For our purposes, therefore, 

F(f,q) = W*)K(q)*-£ (l+| 2 2 ) sin2 ^+^ si n #+.••• (6.28) 

Similar considerations about periodicity and oddness show that 

E(f,q) = P (!-2 2 sin 2 x)Vx=(2/r)E(q)f+J^E m Bin 2mf. (6.29) 

Jo l 

Then, as before, 

(1-g 2 sin 2 ^)^(2/7r)£;((2) + 2X;^^ cos 2mt (6.30) 

i 
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and 



Also 



rir/2 

E m = (2/irm) (1 — q 2 sin 2 x)* cos 2mxdx. (6.31) 

/i 2 • 2 u i <A (2ri-2)!g 2w sin 2n ,. 00 . 

(l- g 2 sm 2 ^=l-g Vfrl(n-l)l * (6 ' 32) 

Inserting (6.32) into (6.31), we find 

^-" ma ^W^Dl J. sin2 "^os2m^. (6.33) 

On inserting (6.24) into this, we obtain 

(-!)»+' ^ (2n)!(2n-2)I^» 

ffm m £j>2 in - l n\(n-l)\(n+m)\(7i-m)\ [ ' 

Thus E tm is of order g 2m or of order k m and 

^ 1= l+i+ ' • •' (6 - 35) 

**--268+ ' • ' ■ (6 ' 36) 

Then, for our purposes 

E(+,q) = (2M E(q) ++(£+£) sin2^^ sin #+.... (6.37) 

Finally, inserting (6.28) and (6.37) into (6.10) and (6.12), we have for the ^-integrals 

Ni and N 2 : 

M = (a 2 -a 3 2 )-^?[5 1 ^i(l+i^sin2^+gsin #+ . . .J, (6.38) 

N 2 =(a 2 -c4)-ho[B 2 +-^ ( 4 + 3 ^ 2 ) sin 2^+|^ sin 4^+ . . .], (6.39) 

with 

5l =2jL 2 [K( i )-E(q)]=l+^+^+ . . ., (6.40) 

B 2 ^-K{q) = \+\ 2 2 +^ 2 *+ . . . . (6.41) 

Here the terms in \p are exact. In N 2 the sines are correct to 0(Jc 2 ), while in Ni they are correct 
only to 0(k); this is as much accuracy as we need for N i} however, since it is multiplied by 
c 2 =kp 2 in the first kinetic equation (2.5). 

6.2. The Integral N 3 
From (6.3) and (3.34) we have 



Then since 
we obtain 
where 



^-^^S^^ (6 - 43) 

L m ^± j^i-^y^i-vyrjiy^drj. (6.45) 
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With 
and 



U=± J\l-v 2 )-\l-V 2 /vir>dv (6.46) 



(l-i? 2 )-V w =(l-i? 2 )- 1 --2> 2B , (w^l), (6.47) 

we then find 

m— 1 

L m =U-^2L ln , (ro^l), (6.48) 

w = 

where 

Li,= ±JV(l-W^)"^. (6.49) 

To evaluate i , rewrite (6.46) as 

^iJV 3 ^- 2 -!)- 1 ^- 2 -^ 2 )"^ (6.50) 

and introduce the new variable X, defined by the equation 

taiix=(l-i;?)* tan^ = |cos/|tan^ (6.51) 

and the requirement that X and \f/ shall keep in step. I.e., whenever ^ equals a multiple of 
tt/2, x shall be equal to \p. Then 

csc 2 i£=l+cos 2 /cot 2 X (6.52) 

and 

rT 2 =77 - 2 csc 2 iA=l+cot 2 /csc 2 X, (6.53) 

so that 

rj~hlrj = cot 2 Icsc 2 X cot X dx (6.54) 

and cot xdx^O according as drj^O in (6.50). With 77 = 7/0 sin \f/ its integrand becomes 

± | tan I\ cot X | tan X | dx= | tan / \dx, (6.55) 

so that 

L =|tani] X=r ?0 (l-^)" i X. (6.56) 

To fit the angle X into one's knowledge of the corresponding elliptic motion, note that as c 
approaches zero, rj approaches sin 0. Then, by (6.53), 

, / tan0 \ ~ ,„ _ wN 

x -^ in \\i^r\) = +- Q > < 6 - 57) 

where 12 is the right ascension of the node. 

To evaluate the integrals Z ln , put 77 = 770 sin \p in (6.49). Then 

L ln=7} 2n+1 r s [tfn X( lx, (6.58) 

so that 

£10=^0^. (6.59) 

To handle the cases n^ 1, we rewrite (6.22) as 

thus finding 

_ ,r'(2nW ^ (-l)'(2n)Isin2# 

jLlre 2 2 *(«!) 2 +"° Z Pi (n+j)\(n-j)\j ' { - l) - {b ' bl) 

Insertion of (6.48), (6.59), and (6.61) into (6.44) then yields 
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{a 2 c4) l\ s -U2-, Q 2 2m (m!) 2 ,0 %& 2 2 "'(m!) 2 ^^ 2 2 »(m!) 2 & 2^(w!) 2 

_ - (2m)!iy 2 - w ^ „2»+i 9 -2„ f, (-l)'(2»)!ain2j> 

If in (6.62) we now use (6.56) and the relation 

(1 "^ } S 2*"(ro!)» ' (6 - 63) 

we find 

(ai-^)W3=W(l-i7§)-*(l-^" 2 )- i X+B3^+S5 8 ,sin2^], (6.64) 



5=1 



where 



and 



£ 3 =l-(l-*7 2 - 2 )-*-;£ 7mv; 2m , (6.65) 

TO = 2 

_ (2m)! ^ (2n)li>g» . 

V^B ,in2,'^- T' ( 2m ) ! ^" 2 " ^' „»»;-..^ (-1) ; '(2tQ! sin 2j> . 

gB„. sin 2, *— g 2 22m(w!)2 g, 2 g (w+i)!(7l _ i)!j (6-67) 



The easiest way to isolate the coefficients B 3s is to use the orthogonality of the functions 
sin 2s\p. Then 

- (2m) l^- 8 - ^ qg"2-»»(-l)'(2n) ! 
^ 3S ~ m f/ +1 2 2 ™(m!) 2 £j (n+«)!(n-«)b ^^ 

= O0h 2s " 2 ), (6.69) 

so that B 3l = 0(v2 i ) = 0(k 2 ), B 32 =0(v2 6 ) = 0(k z ), etc. Thus, to obtain the sine terms of iV" 3 
to 0(k 2 ), we need only 

*«=§ ^ 2 - 4 + .... (6.70) 

oo 

To test the convergence of the series X) Yw*?2~ 2m > note first that 

2 

_gf ! )l_ = 1.3.5...(2n-l) < l 

2 2 *(ra!) 2 2.4-6...2» =2' l = l) ' ^- n) 

and hence that 

r^ (2n)\ v l"+ l m-l 

ti x 2 2 "(n!) 2 = 2 {K)J ~ J 

for any orbit, polar or nonpolar. Then, by (6.66), (6.71), and (6.72) 

Thus r M <m and 

S7m^ 2m <f:m, 2 - 2 -. (6.74) 

ra=2 w=2 

But 

Zj mr? 2 — ,, - 2 \2 ^2 — /-, - 2 \2 ' 1° '-v 



2 



: (l-^ 2 - 2 ) 2 ^ ^ (l-^ 2 ) 2 ' 



where t? 2 2 =0(k). Thus the series S^Ym^ 2 ™ converges rapidly. 
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To test the convergence of the Fourier series Sf=i^3* sm 2s& note that 



But 



'§ ^ Sm 2 **' " § ' Asl " .Si ^rnV S 2«»(n+,)!(n-*)f 



(6.76) 



(2*)! 



(2rQ!yi!n! 



(2n)! (7i)(ri-l)(n-2) . . . (n— *+l) (2n)! 



2 in (n+s)\(n—s)\ 2*"(ra!) s (»+s) !(»—*)! 2 2B (n!) 2 (»+«)(»+*— l)(»+s— 2) . . . (7i+l) = 2 2 "(n!) 2 
by (6.71). Then 



^ l _ (2n)l 

^—-' 92W/ 



and 



t^2 2w (7i+^)!(7i-5)! 
S^3 S sin2^ 



<^m 



< £ S- ta <s **- 



w=s+l 



2 2m (m\) 2 



by (6.71) again. Then, by (6.75) 



S A. sin 2,0 






7? 2 4 (2 — r? 2 2 ) 



(l-* 1 ) 1 



(6.77) 



Since r)2 2 =0(k), the Fourier series converges rapidly. 

Thus all terms in (6.64) remain finite for all r] 2 ^ 1, except the term involving x, which appar- 
ently may become infinite for polar orbits. Note, however, that N 3 occurs only in eq (2.7) 
for the right ascension, when it lias a factor a 3 , which vanishes for a polar orbit. To see what 
happens in this case we must investigate the limit of a 3 N 3 as %~^1. For a polar orbit we have 

a z N 3 =a 3 (c4-la)^ Vo (l-vl)-K^-V2 2 )-^+0(a z ). (6.78) 

But \a- s \(c4-al)->=\a 3 \a; 1 to 0(a 3 ) and by (4.8) and (4.2) 

kk 2 - 1 =(l-i7?)Kl-^ 2 sin 2 i^=(l-v^) h (l-V2 2 ) 1 to 0(* z ). (6.79) 

Tims to 0{a 3 ) 

a z N z =(sgn a 3 )x (polar orbit). (6.80) 

If we now use (6.51) to plot x versus ^for various values of 775, we obtain figure 1. 




Figtjke 1. Plot of x versus \p for various values of 



7o« 
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The general form of these curves can be checked, as follows. From (6.51) 

tan X=e tan & e = (l — ij§)i (6.81) 

To find what happens to a polar orbit we must let € approach zero. 
If we put 

*=(n+J)ir+Aifc X= (»+*)*■+ AX (w=0,l,2.3, . . .), 

then 

tan Ax=€ -1 tan A^. 



7T 



If we keep A^ fixed and let € -»0, then if 0<AiK>/2, it follows that A x ->7r/2. If -£<A^<0, 

then Ax— >— 7r/2. In a polar orbit, whenever the satellite passes over a pole, ^=(^+i)?r, 
with ^>0? so that % thus jumps by +ir. By (6.80) and (2.7) the right ascension <t> then jumps 
by +7T in a direct orbit or by — t in a retrograde orbit. These are expected results, which 
had to be obtained as a partial check of eq (6.64) for N z . 

7. Mean Motions 

The purposes of this section are to obtain expressions for the mean motions correct through 
the first order, to find if we are on the right track, and to obtain exact expressions for the mean 
frequencies for later use in section 8 in checking the secular terms in the final solution. 

If p p , pr,, and p<p are the generalized momenta conjugate to p, rj, and <£, then the action 
variables 

Ji=(yp4p=2jPidp, (7.1) 

J 2 =(£p r] d7 1 = 4\" p v dr 1 , (7.2) 

Jz=(vpd4>=\ p<j>dc/>=2ira 3 , (7.3) 

are functions of the Jacobi constants a u a 2 , <* 3 . (Since these e/ ? s occur only in this section, 
there is no danger of confusion with the coefficients of the zonal harmonics in the expansion 
of the potential.) The a's are then functions of these J's and the mean frequencies [8] are 
given by 

mean p-frequency = v x = daJdJ u (7.4) 

mean r7-frequency=v 2 = da:i/dJ2, (7-5) 

mean ^-frequency = *> 3 = daJc)J z . (7.6) 

Note that v 2 and *> 3 are identical with the usual nodal and sidereal frequencies, but that v x is 
somewhat different from the usual anomalistic frequency. 

To compute these frequencies, one may use the system of equations 

23 xr~ 7s^- = ?^r =8ln > (w= 1,2,3). (7.7) 

m=i OJ m Oa n Oa n 

With the use of (7.4) through (7.6) and the abbreviation 

Jmn = Z>Jmlda n , (7.8) 
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these equations become 

^1^11 + ^21 =1, (7.9) 

viJi*+*Jn =0, (7.10) 

^ie7i3 + ^2e/23 + 27n/ 3 =0, (7.11) 

with the solution 

^=J 22 /A (7.12) 

v 2 =-J 12 /A (7.13) 

27tv s = — viJ iZ — V2J2Z (7-14) 

A = JnJ22 — 1/12^21 v'.loj 

With the use of eqs (13), (53), and (55) of [1] and sections (5) and (6) of the present paper 
we find 

J 11 =2Si(p 2 )=27r(~2a 1 )-*(a+6i+^i), (7.16) 

J 12 = _2a 2 /i > 2 (p 2 )--27ra 2 (-2a 1 )-^ 2 , (7.17) 

J n =2cZath\(p. 2 ) = 2KC*a 3 (-2a l )-*A 3 , (7.18) 

J 21 -4c 2 M(W=27rC 2 (^-al)-^2/^ (7.19) 

J 22 = 4a 2 N 2 (r l0 )=2Ta 2 (al-ai)-> V{) Ii 2 , (7.20) 

J23= -4a z N^ Vo ) = -2ira 3 (^-^)-iT ?0 [^3+(l-i?2)- i (l-'?2" 2 )"*]. (7-21) 
Then 

27r^i=(~2a 1 )i[a+6i+A+c 2 ^ 2 5 1 B 2 " 1 ]- 1 (7.22) 

27r? 2 =(a 2 -a 3 2 )*7 7o - 1 ^ 2 ^ (7.23) 

The above results are all exact, for the potential (2.3), and we shall refer to them again in 
section 8. 

For comparison with other theories it is desirable to express these mean frequencies as 

power series in the oblateness parameter. For this purpose it is convenient to use the orbital 

elements a 0) e , and i , and since we shall carry the series only through the first power we may 
replace k = c 2 /pl by k = c 2 /p 2 . 

With use of the relations in sections 3 and 4, the above equations then become 

2irv 1 =n +O(k 2 ) (7.24) 

2wv 2 =n [l + ik(5 cos 2 i -l)] + O(P), (7.25) 

H=nlal (7.26) 

J lz =irk(2+el) cos i Q + 0(k 2 ), (7.27) 

J 23 =-27rsgn a,+wk(l-el) cos i +0(k 2 ), (7.28) 

which together with (7.14), (7.24), and (7.25) lead to 

27rv 3 K=(sgn <*)[l+*fc(5 cos 2 ?: -l)]-f k cos i +0(k 2 ). (7.29) 

Here sgn a 3 =±l accordingly as the orbit is direct or retrograde, respectively. 

To avoid any use of the concept of an osculating ellipse, we may define the mean motions 
as follows. We say that the ascending node exists only when the satellite is over the equa- 
tor, travelling north. Let 12 be its right ascension at such a time. We then define the mean 
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where n is given by 
Similarly 



motion of the node relative to OX by 

5-lim 5i=5?, (7.30) 

where Q and £ are the values of ft and t at some (ascending) node and ti t and £< are their values 
i nodes later. Since the present system is of the conditionally periodic Staeckel type, it 
follows, after some fairly close reasoning that I shall here omit, that 

5=ar(|Fa|-i>a) sgn <*. (7.31) 

We also say that p-perigee exists when, and only when, p=pi and we let <£ be its right as- 
cension at such a time. Then the mean motion of p-perigee relative to OX must be equal to 

the mean motion <£ of its equatorial projection relative to OX and 

I=lim fc?- (7.32) 

Here $o and t are the values of <i> and t at some p-perigee and $* and t% are their values i 
p-perigees later. Then, again omitting the proof, we have 

i=27r(|^|-^i)sgna 3? (7.33) 

Now let co be the arc on the celestial sphere from an ascending node (when it exists) to 
the next p-perigee. Then, since the mean relative motion of these points must be equal to 
the mean relative motion of their projections in the equatorial plane, we find that the mean 

motion co of p-perigee relative to the node is given by 

Z=(l-V)sgna z , (7.34) 

=2t(v 2 -p 1 ) ) (7.35) 

with use of (7.31) and (7.33), for any orbit, direct or retrograde. 
From (7.31), (7.25), and (7.29) it follows that 

tl=-ikn cos i +O(k 2 ) (7.36) 

and from (7.35), (7.24), and (7.25) that 

l=ikn (5 cos 2 i -l) + O(k 2 ) (7.37) 

Here n /2ir is the frequency in an elliptic orbit with the same total energy. Equation (7.36) 

agrees with results found by many other authors, as does (7.37) when co is the mean motion 
of r-perigee relative to the node. By (2.1) and (2.2), however, r 2 = p 2J rkp 2 (l — rj 2 ), so that r 
and p differ by a variable quantity of 0(k). It is thus a little surprising that the mean mo- 
tions of p-perigee and of r-perigee relative to the node should be equal through 0{lc). This 
means that the mean p-frequency is equal to the mean r-frequency to this order. 

8. Solution of the Kinetic Equations 

Before solving the kinetic equations (2.5) through (2.7) it is convenient to have several 
relations connecting the uniformising variables E and v. From (5.12) we obtain 
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cos E—e . v 

cos#=- ^ (8.1a) 

1 — e cos E v 

The requirements that dv/d£>0, dE/dt^O for all £ lead to the result that dv/dE^X) for all t. 
Because of this result, (8.1a) leads to 

sin^+^j =r> (8.1b) 

l — e cos £; 

without ambiguity in sign. For a given value of E } eqs (8.1a,b) determine v modulo 2-k. 
On imposing the further requirement that v shall always equal E whenever the latter is a 
multiple of ir, we find that E determines v completely. Two other relations arc often useful, 
viz, 



tan ^ 



and 



=(S)* tan f (8 - lc) 



v—E 7 sin E _ ir f oxn/, /o-u\ 

^X Vtcob^ y " e [1_(1 "' )jI<1> (S - ld) 



Before beginning the solution of the kinetic equations it is desirable to assemble the 
results already obtained. By (2.5), (2.6), (5.1), (5.2), (5.30) through (5.33), (5.6), (5.28), 
(5.35) through (5.40), (6.1), (6.2), (6.6), and (6.38) through (6.41), 
the equations for p and 77 are 

t+p l =(-2a l )- 1 ^[b 1 E+a(E-e sin E)+AiV+A u sin v+A u sin 2v] 

+c\a\-c$)-Hl ["-Bi*-4 (2 + g 2 ) sin W+~ sin #l+periodic terms of 0(F), (8.2) 

&/a&= — (— 2ai)~?[A 2 v+A 2 i sin v-\-A 22 sin 2v+A 23 sin 3v+A 24 sin 4v] 

+ (at-ai)~^ [b^-^ (4+3<? 2 ) sin ty+^jfc sin 4^1+periodic terms of 0(F). (8.3) 

Here 

p = a(l— e cos E)=a(l— e 2 )(l+e cos i?)" 1 , (8.4) 

V=Vo sin ^, (8.5) 

with E and p connected by any of the eqs (8.1). 

A 1 =(l-e^p ± (6 2 /^) w P,(6i/62)^,-2[(l-^ 2 ) i ]=0(«, (8.6) 

71 = 2 

A =(l-e 2 )^ 1 i; (6 2 /p)"P K (V6 2 )BJ(l-e 2 )5]=0(fc ), (8.7) 

A ll =Kl-e 2 )ip-h(~2b l b!p+bi) = 0(k 2 ), (8.8) 

A 2 =^ (l-e*)»2T 8 &2e*=0(fc*), (8.9) 

fli=^[*<fl)-.E(2)]=0(*"), (8.10) 

5 2 =?K( 2 ) = 0(F), (8.11) 

a=flo/V2=0(ki). (8.12) 

-fiT(<Z) and i?(g) are the complete elliptic integrals of the first and second kinds. 
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A 21 =(l-e*)hp- 5 \-\Af+(^A*-By p *+\ ABp (i+0+| J 52(4+3€ 2 )]-O(/:), (8.13) 
A 22 =(l-e>)*p-* [(|^ 2 -^) ~P 2 +^ ABe>p+l B 2 Qe 2 +\e^y=0(k), (8.14) 

A 23 =(l-e^p- 5 e £^ + ^=0(*»), (8.15) 

^ l21= 2§6 ^ l -^) h V- 5 B^=0{¥). (8.16) 

In the above equations c 2 =r 2 e J 2 is regarded as known, as are the orbital elements a, e, 
7)o, fa, and ft,. Then p=a{l-e 2 ) and rj 2 , A, B, on, a 2 , and a 3 are given by section 4. Also 

&i=-|a, 61=5, R m (x)=x m P m (\lx), with |aj|<l. 
To solve (8.2) and (8.3), place 

E=E S +E P , v=v 8 +v p , t=t s +t P , (8.17) 

where the subscript s means "secular" and the subscript p "periodic." Then if p goes through 
Ni cycles in time T x and if 77 goes through iV 2 cycles in time T 2 , 

£ = l = E s==is== i im ^=2^1 (8.18) 

£=^=lim 2 -^=2ttv 2 . (8.19) 

Since we have already obtained exact expressions in section 7 for v x and v 2) it is clear that we 
can obtain the secular terms exactly. We shall also obtain the periodic terms through 0(k 2 ). 
By (8.18) we can write 

E s =v s =M s , (8.20) 

where M s will play the role of the secular part of a mean anomaly. 
Then 

E=M,+E P , v=M s +v p , *=*+*,- (8.21) 

We may obtain the secular solution of (8.2) and (8.3) independently of section 7, by 
dropping all the sines in those equations, placing E=v=M s and yp=^ s , and solving the resulting 
equations for M s and \p s . The resulting equations are 

(-2a 1 )-Ha+ h+AJMs+c 2 (a 2 -al)-hlB^ s =t + ^ (8.22) 

-(-2 ttl )-* A 2 M s +(a 2 -ci)-hoB 2 & = P 2 la 2 , (8.23) 

with the following solution. 

8.1. Secular Solution 

M ( 9 V B 2 (t+fi l )-C 2 p 2 a 2 ^ 2 B 1 , 

(a+bi+Ai) B 2 +c 2 t)lA 2 B x 

, ,2 2 \* -1 A.jt+Pd+faa^ia+^+Ar) 

+ s =(a 2 -a s y Vo (a+bl+Al) B 2+ ^ lAaBl ' (8-25) 

Comparison of these results with (7.22) and (7.23) shows that ilf s =27n< l and \p s =2irv 2 , as ex- 
pected. We may now rewrite these expressions more conveniently as 

M s =2Tv 1 (t+fi l -c^ 2 ^ l v lB l B^), (8.26) 
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t s =2irv 2 [t J rPi J rP2^ 1 {a+b 1 +A 1 )A^}. (8.27) 

As a check, note that to 0(k°), 2irv l = 2Trv 2 =n !t) so that 

2lf.=»io(*+A)+0(fc) l 

^ =n ° (HA)+ S +0(t) ' 

But, to order k°, a 2 =(np)? and A 2 = {l-e 2 )*jr l , so that ntfi(a 2 A 2 )~ 1 =n a2)[n2)(l-e 2 )]~l=l, 

since nla?=ix to 0(k°). 

Thus 

h=M s +&+0(k), 

as is to be expected, with /3 2 replacing co. 

As a later aid in reducing the solution for the periodic terms to Kepler's equation, it is 
here convenient to rewrite (8.22), by transposition of some of its terms: 

t+p 1 -(-2a 1 )~Ka+b 1 )M s =(-2a 1 )-*A 1 M 8 +C 2 ^ (8.28) 

8.2. Periodic Terms 



We shall put, successively, 








E p =E 0f 


V p =Vo, 


V^=1A>; 


(Step 0) 


E p =E -\-Ei, 


V p =Vq+Vi, 


^p=^0+^U 


(Step 1) 


E P =E Q +Ei+E 2 , 


Vp=V Q +Vi+V 2 , 


^p = lfo+lh + lfe. 


(Step 2) 



In step 0, we retain in the equations (8.2) and (8.3) only those periodic terms which are 
of (Ok ), viz, sin E. In step 1, we retain in these equations all periodic terms of order k° or k, 
but none of higher order. In step 2, we retain in the equations all periodic terms of order k°, 
k, or k 2 , but none of higher order. In carrying out each step, however, we shall suppose that 
each quantity involved is calculated to such an accuracy that the error is of order k 3 . Then, 
effectively, E , v , and \p will all contain terms of order k and k 2 , as well as terms of order k°. 
E u v 1} and \pi will contain no terms of order k°, but will contain terms of order k and k 2 . E 2 , 
v 2 , and \p 2 will be of order k 2 . Such a procedure will greatly simplify the resulting equations. 

8.3. The Periodic Contributions E , v () , and \p Q 

On placing E=M S +E , v=M s +v , and ^=^ s +^ in (8.2) and (8.3) and retaining only 
the term sin E of the periodic terms, we find 

t+fi 1 =(-2a 1 )^[(a+b 1 )(M s +E )-ae sin (M s +E )+A 1 M s ]+c 2 (c4-a 2 )-^lB^ 8) (8.29) 

ft/^2=-(--2aO-^4 2 (M s + ^o) + (^-a3 2 )-^o5 2 (^+^o). (8.30) 

On subtracting (8.28) from (8.29) and dividing the resulting equation by (a+bi)( — 2a 1 )~*, 
we find 

M,+E -e' sin (M,+E Q )=M 8 , (8.31) 

where 

ae 



a • />, 



<1, (8.32) 



since &i>0. Equation (8.31) is Kepler's equation for M s +E , with an effective eccentricity 
e'. Let us suppose it to be solved by the most appropriate method, which will depend on the 
value of e'. We then have M s -\-E and can then find v=M s J r v by use of eqs (8.1). 
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On subtracting (8.23) from (8.30) we then obtain \f/ Q as a function of v$\ 

^={-2a l )-Kc? 2 -c$)ho l A 2 B;W (8.33) 

Here the coefficient \p fv is unity to 0(k°), but to follow the procedure outlined above we must 
not make such an approximation. Instead we must calculate it so accurately that the error is 
of order k z . 

8.4. The Periodic Contributions E u vi 9 and \pi 

Now, knowing M 8 , \p s , E , v , and \p Q , we place E=M s -\-Eq-\-Ei, v=M s -\-v Q -\-v u and \p= 
^As+^o+^i into eqs (8.2) and (8.3), discarding only the periodic terms of order k 2 . Then 

t+^=(-2a0-'i(a+b0(M s +E o +E 1 )-ae^m(M s +E o +E l )+A 1 (M s +v o )] 

+ c 2 (at-a!)-hl[B 1 (rP s +^)-ism (2&+2* )], (8.34) 

2 /a 2 =-(-2a 1 )-*[A 2 (M s +v o +v 1 )+A 21 sin (M 8 +v ) + A 22 sin (2M s +2v )] 



+ («S-ai)-*i ? o[-B 2 (^+^+^)-f sin (2fc+2* )]. (8.35) 



Subtraction of (8.28) from (8.34) and division of the result by (a+b i )( — 2a l )-? now gives a 
Kepler equation for M s +E +Eii 

Ms+Eo+Ei-e' sin (M S +E +E 1 )=M S +M U (8.36) 

where 

Afi-Ca+W" 1 [ -(A 1 +c^lA 2 B 1 B 2 l )v Q +j (-2^)^1-4)- hi sin (2*,+2* )]. (8.37) 

Here we have used (8.33) to combine terms in v and i/v The quantity M s -\-M\ is then a mean 
anomaly whose secular part is exact and whose periodic part is correct through order k. It 
has no periodic part of order k°; this is characteristic of a mean anomaly. 

It is not necessary to solve the Kepler equation all over again. If in (8.36) we put 

sin (M s +E +E 1 ) = (l-±Ei) sin (M s +E )+E 1 cos (M s +E ) + O(El), (8.38) 

the error is of order k z . Then (8.36) and (8.38) yield a quadratic equation for E u whose so- 
lution through terms of 0(k 2 ) is given by 

F _ Mi e' M\ sin (M s +E ) 

1 \-e' cos (M s +E ) 2 [\-e f cos (M s +E )f ' } 

To find Vi insert v=M 8 +v +Vi and E=M S +E +Ei into eqs (8.1) and solve for v x . 

On subtracting (8.23) from (8.35) and eliminating terms in v and \j/ by use of (8.33), we 
then find 

^ 1 = (-2a 1 )--(a 2 2 -a 3 2 )^ ~ 1 -B 2 - 1 [^i+^2i sin (M s +v )+A 22 sin (2M s +2^ )] (8.40) 

+| 2 5 2 - 1 sin(2^+2^ ). 

(Note that the elimination of v and \[/ Q would not have been possible if they had been carried only 
through order k°.) 

8.5. The Periodic Contributions E 2 , v 2 , and \[/ 2 

Finally, knowing M s , \// s , E , v , \p Q , E ly v h and \pi, we place E=M s -\-E -\-Ei J rE 2 , v=M s 
+^0+^1 + ^2, and ^^tAs+^o+^i + ^2 in (8.2) and (8.3), discarding only the periodic terms of 
order greater than k 2 . The equations become 
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t+l3 i = (~2a l )-i[(a+b l )(M s -\-E +E 1 +E 2 )~aesm(M s +E +E 1 +E 2 )+A 1 (M s +v +v 1 ) 
+A n sm (M,+v )+A 12 sin (2M s +2v )]+e i (ai-al)-i v s jB l (^+^+4' 1 ) 

-~ (2 + 2 2 ) sin (2*,+2*«+2* 1 )+|£ sir. (#.,+#„)], (8.41) 

/3 2 /a 2 =-(-2a 1 )-*[^ 2 (M J +»o+«i+»2)+^2i sin (M s +o +»i)+^22 sin (2M s +2p +2v l ) 
+A 23 sin (3M,+3» )+^24 sin (4M s +4v Q )} + (a 2 2 -c$)-i Vo [ J B 2 (*+*„+* I +* 2 ) 

-g (4+3«f) sin (2^+2*0+2*0 +Jg <Z 4 sin (4&+4* )]- (8.42) 

On subtracting (8.34) from (8.41) and discarding periodic terms of order k 3 or higher, we find 
0=(-2a 1 )-i[(a+b 1 )E 2 -aeE 2 cos (M,+E +E 1 )+A 1 v 1 +An sin (M s +v ) 

+A 12 sin (2M s +2* )]+c 2 (^-a|)-^ [Wi-^ h cos (2*,+2* ) 

(2*+2*,)+£ sin (4*,+4*o)l- (8.43) 



2 2 
— tt sin 



Then 



where 



Ei 1-e' cos (mI+Eo+EJ (8 - 44) 



M 2 = - (a+Si)- 1 ^!?;^^! sin (M,+fl ) +A 12 sin (2M s +2t>„) 

+c i (-2a l ) i ( a l-c4y i r,lf B^-^h cos(2*+2* )-| 2 sin(2*+2*„)+£sin (4&+4* ) !]• (8.45) 

It is easy to show that, to order k 2 

Ms+Eo+Et+E.-e'sm (M s +E +E l +E 2 )=M S +M 1 +M 2 , (8.46) 

so that M 2 is the second-order periodic term of a total mean anomaly 

M=M S +M 1 +M 2 + . . . (8.47) 

corresponding to the effective eccentricity e'. 

To find v 2 insert v=M g +v +v 1 -\-v 2 and E=M,+E +Ei+Ei into eqs (8.1) and solve 
for v 2 . 

To find \f/ 2 , subtract (8.35) from (8.42) and discard periodic terms of order k 3 or higher. 
The result is 

* 2 =(-2ai)-Ha|-«3 2 )^o" 1 ^2"U^2t' 2 +^2i%cos(M s +»o)+2^ 2 ^iCOs(2M s +2y )+^ 2 3sin(3M s + 3^) 
+^ 24 sin(4M s +4^)]+-J J B 2 - 1 ^ 1 cos(2^+2* )+3| 2 sin(2* 2 +2* )-3£sin(4*+4* () )J- (8.48) 

This completes the solution, with exact secular terms and periodic terms correct through 
order k 2 , for E, v, and * and thus for the spheroidal coordinates p and -q. 
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8.6. The Right Ascension <£ 
From (2.7), (5.3), and (6.3) we obtain 

<t> = fc + aM-c 2 *zRz. (8.49) 

Then, from the equations in sections 5 and 6 

0-ft+a 3 (^-a3 2 )-^ o [(l-^)-Kl-r; 2 - 2 )^X+53^ 

+;§ ^2~ 4 sin 2*1 -c 2 a 3 (-2ai)-^ [^ + S ^a» sin wl, (8.50) 

where x is an angle that always equals * whenever * is a multiple of 7r/2 and which also satisfies 

tan x=(l — rjl) 1 * tan *. 

The expressions for A z and the Az n 's are given in (5.61) to (5.65) and that for B z is given in 
(6.65). With the secular parts of v and \p exact and their periodic terms correct through order 
k 2 , the right ascension 0, as given by (8.50), has a secular part that is exact and a periodic part 
correct through order k 2 . To check the secular part of </>, note that one can obtain it from 
(8.50) by placing x=\f/=\p s and v=v s and discarding the sines. If we do so and also use b s =2wvi 
and \[/ 8 =2ttv2, we find 

* f =-2*ft*(-2<^-M«iri+^^ (8.51) 

on comparison with (7.18) and (7.21). Thus we find 

J>=4> s =2ttv 3 (8.52; 

by (8.51) and (7.14), a result known to be correct [8]. 

A summary of the principal results of the paper follows in section 9. 

9. Summary of the Solution 

We assume that /jl and c are known, where n is the product of the gravitational constant and 
the mass of the planet and where c 2 =r 2 e J 2 , r e being the equatorial radius and J 2 the coefficient 
of the second zonal harmonic of the planet's gravitational potential. For the earth e/ 2 ~ 
(1.08)10- 3 . 

If X, Y, Z are the usual rectangular coordinates of an artificial satellite and if r, 6, 4> are 
respectively its planetocentric distance, declination, and right ascension, its oblate spheroidal 
coordinates p, r?, <t> are given by 

X J r%Y=r cos 6 exp % 0=[(p 2 ^-c 2 )(l — rj 2 )]* exp i4>, 

Z=r sin d=pr), (— lg^gl). 
The potential 

F a =- M p(p 2 + cV)~ 1 

then fits the even zonal harmonics exactly through the second and, in the case of the earth, 
approximately through the fourth. Solution for the motion with such a potential thus furnishes 
a very accurate intermediary orbit. Since this potential leads to separability of the Hamilton- 
Jacobi equation, the solution is given implicitly by the quadratures of eqs (2.5) through (2.7). 
The integration constants are the Jacobi as and /3's. 

If the initial conditions are known, one can readily evaluate the as. Then if one can 
evaluate the integrals in (2.5) through (2.7) one can also evaluate the /3's. Evaluating the 
integrals depends on factoring the quartics F{p) and 67(77). The factoring of G(ri) is immediate, 

196 



since it is quadratic in r? 2 . To discuss the factoring of F(p) we introduce, in place of the as, 
the orbital elements a = — p/2a h e = [l+2aial/fx 2 ]^ and i =cos~ l (a 3 /a 2 ). 
With p varying in the range pi^ p^ p 2 , we write 

F(p)^c^l+(p 2 +c 2 )(-al+2pp+2aip 2 ) = (-2ai)(p-pi)(p 2 -p)(p 2 +Ap+B) 

and find that A, B, a, p, and e occur in the p-integrals. 
Here 

a=U Pl + P2 ), i- e .- ? ^ ? ^ a(1 _ e 2 )= ^L, ,.y=£!. [1 _ (1 _^ ] » 

^ (Pl + p 2 j" P1 + P2 P2 + P1 

Thus ^4, 7?, Pi + P2 ? and pjp 2 are needed. If we put 

£>o=&o(l — £0), h=C 2 /Po, x = (l — effi, y = a^a 2 =^cosi , 

then through 0{k 2 ) 

A=-2k p G if[l+h(2x*-3xY-±+Sy 2 )+ • • .], 

B=k oP( ?(l-if)[l+h(4if-xy)+ . . .], 

2a^p 1 + p 2 =2^ x- 2 [l-^y-^Y(2^ 2 -3xY-4+8?/ 2 )+ . . J, 

a^=pip 2 =jir- 2 [l+^7 2 (^ 2 -4)-% 2 (12x 2 -x 4 -20xV-16+32?/+xV)+ . . .], 

If a , e , i , ft, ft, ft are used as orbital elements, we then assume that the quartic F(p) is fac- 
tored numerically, to as great an accuracy as may be desired, with the aid of the above relations 
as a starting point. 

With rj varying in the range— 1^ — rj ^ 7?^ t? () ^ 1, we write 6^77) =— a|+(l— 7? 2 )(a 2 ^2aicV) 
--2a 1 c 2 () 7 2 -7 7 2 )(7 7 2 -r 7 2 ) and find 



WV 2(a 2 -a 3 2 ) L \ (a 2 -2a^) 2 / _[ 



Here t? 2 :$>1. Then all the quantities a, p, e, A, B, rj , and r} 2 are known in terms of the orbital 
dements a , e , and i . 

If we assume that the orbital elements are to be determined by an iterated least-square 
fitting of the solution to many revolutions in the orbit and not by initial conditions, there is a 
better set of elements, introduced by Izsak [4], These, viz, a, e, i^sin -1 ^ ft, ft, and ft, al- 
though not easily found from the initial conditions, result in immediate factoring of the quartics 
F(p) and G(rj). We therefore give the solution in terms of these quantities, with the under- 
standing that they are to be determined either by the least-square fitting or from initial con- 
ditions by numerical factoring of F(p). 

Given p, c, and the elements a, e, I, ft, ft, ft, 
compute 

?7o=sin/, 

p=a(l-e*) } D=(ap-c 2 ) (ap-c'rjD+U^rjl D'=D+4a?c 2 (l-ril), 

A=-2ac*D- 1 (l-r 1 l) (ap-c^t), B=cW,T)~ l D f , 

b, = -\A, b 2 =B$, -2ai=p(a+bi)-\ 

~lr-=<loPo=--C>(l-r,l)+apD- l D', a 2 =(-2ai)KdoPo)<>0, 

a 3 =a 2 ^l-— J COS /, V2 =^#7' k=C i /p 2 , q=Vo/ri2- 
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Restrict considerations to the case b x jb 2 1k 1. Then 

where I c =l°54:', approximately. Equatorial and almost equatorial orbits are thus ruled out. 
Then compute 

A l =(l-e*)ip± l (b i /p)»P tt (b 1 /b 2 )R n - 2 [(l-e*n 

71=2 

oo 

A^il-e^p-^ihJpyP^/bJRAil-ey), 

71 = 

where P n (x) is the Legendre polynomial of degree n and where B n (x)=x n P n (x~ l ), always a 
polynomial of degree [n/2] in x 2 . 

A3=(l-e*)*p-*±,D m R m+2 [(l-e*)i], 

m=0 

where 

A«=S(-l) < -"(c/p)»'-»-(6 JI /p)*»P IJi (6^ 9 ) J 

71 = 

D 2i+1 =±(-iy- n (c/ P r-^(b 2 /pr+ i p 2n+1 (b J /b 2 ) ) 

71 = 

B 1 =2r- 1 q-HK(i)-E(q)] ) B 2 =2^K(q), 

where K(q) and E(q) are the complete elliptic integrals of the first and second kinds,[respectively. 

R i /, -,m ^ -» u (2m)'- <£? (2n)\ v l n 

Ba=l-(l-^ 2 )»-2T^-, where IW-j^pig^ 

(The above series all converge rapidly.) 

An=^(l-e*)*p-'e(-2biHp+bS), Au=^p-(l-*)*bU*, 

il M =(l-e»)»j»-> [I (3&»-6»)p-»-| eikftfr-i+JL 6 2 (6e 2 +e 4 );?r 4 ], 

^4=2§6 (l-« , )*2»- , *J« 4 . 

^ 3 i=(l-e 2 )^- 3 e [2 + 6^-' (3+fe 2 )-^- 2 (J^+c 2 ) (4+3« 2 )], 

^=(l-e 2 )^- 3 [£+f bip-V-pr* (j+| 2 ) (| J i+ c2 )} 
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2irv l ={-2a l )Ka + b l +A l +cSlA 2 B l B 2 ')~\ 

2irv 2 =(a 2 2 -at)ho 1 A 2 B 2 -\a+b 1 +A 1 + c'r } 2 A 2 B 1 B^)~K 

The uniformising variables E, v, and yp are then given by E=M s -{-E p , v=M s J rV p , and 
\l/=ip s J r\[/ p . If t is the time, their secular parts M s and \l/ s are given exactly by 

M 8 =2irp 1 (t+Pi-c%a 2 - 1 r, 2 Q B 1 B 2 1 ), 

rP,,=2TV 2 [t+^+^a 2 1 (a+b 1 +A 1 )A 2 1 ]. 

Let the periodic parts be split as follows: E P =E +Ei+E 2j v p = v -\-v ] + v 2 , and 
\l/ p =ypQ J r \pi J ryp 2} where, e.g., E contains terms of order k°, k, and k 2 , E Y contains terms of order 
k and k 2 , and E 2 contains only terms of order k 2 . 

Then E is given by the Kepler equation 

M a +Eo-e' sin (M S +E ) = M S , 

where e f =a(ae J rb l )~ l <^e. The term v is then given by placing ^=il/ s +z? and E=M S +E in 
the anomaly connections 

cos ?;=(cos E—e) (\ — e cos E)~ l sin #= + (1 — e 2 )*(l — 6 cos E)~ l sin fi 

or equivalent relations. (Note that the e here is the original e and not the e' in the Kepler 
equation.) Then 

*>= (-2*)-* (oj-oj)* T70" 1 4 2 B, 1 Vq. 

The term ^ is now given by 

E 1 =[l-^cos(M s +^)]- 1 M 1 -^ / [l-6 , cos(M s +^ )]" 3 M?sm(Af s +£ , ), 



where 



Mi = (a+bJ- 1 [-(A 1 +*rio*A 2 B l B 2 - l )v +^ 



The term t? x is then given by placing v=M s +v +Vi and E=M s -\-E -\-Ei in the anomaly connec- 
tions. Then 

^ 1 ==(-2 ai )-*(ai-^)^o" 1 5 2 -M^i+^2isin(M 8 +!?o)+A 22 sin (2M s +2^ )]+|^ 2 - 1 sin(2&+2* ) 

Finally 

# 2 =[l-6 r ^(M.+^o+Si)]" 1 ^, 
where 

M^-ia+bJ-'lA^+An sin (M,+i?„)+A2 sin (2M s +2%) 

+c 2 (-2*)* W-Q»-*iyi! / Aft-~ A cos (2&+2&)-|sin (2&+2^ )+£ sin (4&+4* ) \]. 

Then ?; 2 is found by placing v=M s J rv -{-Vi-\-v 2 and E=M s -\-E -\-E 1 ^-E 2 in the anomaly con- 
nections and 

^ 2 =(-2a 1 )-K^-^ 2 )^ - 1 J5 2 - 1 [A 2 ^ 2 +4 21 ^ cos (M s +z> ) 

+2A 22 V! cos (2M s +2%)+^ 23 sin (3M,+3tf )+4 24 sin (4M s +4^ )] 

+^ 5 2 - J [^ cos (2^+2^ )+^- 2 sin (2&+2^ )~^ sin (4fc+4* )]- 

The spheroidal coordinates p and 77 are then given by 

p=a(l — e cos J?) = (1+6 cos v)~ 1 p, v = Vo sin ^, where 
E=M.+E +E 1 +E 2 , v=M s +v +v 1 +v 2 , *=*. + * +*i + *2. 
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The right ascension is 

4=h+<*(<4-ci)-*i,o [(l-^)-Kl-^" 2 )- 1/2 X+5 3 ^+^ rjlrj^ sin 2*] 

— c 2 a 3 (— 2qji) _ * A 3 «+2 ^3n sin nv I 

Here X is an angle that equals \p whenever \p is a multiple of 7r/2 and which also satisfies 
tan X= (l — rjo)^ tan \f/. If a? 3 is positive or negative, the orbit is respectively direct or retrograde. 
The above solution gives secular terms of the intermediary orbit exactly and periodic 
terms correctly through order k 2 . There are no long-period terms. 



10. List of Symbols 

We here list only the symbols that are most frequently used, giving for each a short defini- 
tion or the number of the equation in which it first appears. Note that the first digit in each 
such number is the number of the section. 

Definition or equation number 
-action variables; (7.1) through (7.3), 
—7 8 

-3.16; cVpl 
-4.17; c 2 //> 2 , 
-6.15; complete elliptic integral of first kind, 

modulus q, 
-6.46, 
-6.45, 
-6.49, 

-8.20; secular part of "mean anomaly," 
—8.37, 8.45; periodic parts of "mean anom- 
aly," of orders k and k 2 , 
-6.1 through 6.3; 77-integrals, 
-aoil-eo), 
-a(l-e 2 ), 

— Legendre polynomial of degree n, 
-1/0/172, 

—2.1 and 2.2; geocentric distance, 
-equatorial radius, 
-3.1, 3.2; zeros of f(p), 
—5.1 through 5.3; p-integrals, 
-r*PJjr l ), 
-sign of a 3 , 
—time, 

—5.12; "true anomaly," 
—8.17; secular part of v( = M s ) 
—8.17; periodic part of v, 
—various terms of v v , 
—2.3; the potential of this paper, 
—2.1 and 2.2; rectangular coordinates, 
-3.17; (1-^) 1/2 , 
—3.18; a 3 /a 2 = COS to, 
—total energy, first Jacobi const, <0 for a 

satellite orbit, 
—second Jacobi constant, 
— Z-component of angular momentum, 

third Jacobi constant, 
—Jacobi /3's, 
-8.1 d, 
-6.66, 
-5.47, 
-6.81, 



Symbol Definition or equation number 




Symbol 


a — 3.23; element of 2d set, 




J 1, J 2, J 3~ 


a — 3.3; element of 1st set, 




J mn~ 


4—3.11 (also 4.12), 




k - 


Ayr- 5.31, 




k- 


A 2 — 5.36, 




K(q)- 


A u — 5.32, 






A n — 5.33, 




Lr 


A 2] — 5.37 through 5.40 for j=l, 2, 3, 4, 




L m - 


4,— 5.61, 




L\n 


A zi — 5.62 through 5.65 for j=l, 2, 3, 4, 




M 8 - 


biM — 5.6, 




M h M 2 - 


B— 3.11 (also 4.13) 






A— 6.40, 




N u N h N r 


#2— 6.41, 




Po- 


£ 3 — 6.65, 




P- 


£ 3 — 6.64 (also 6.68) 




Pn 


c — 2.4 (fundamental distance in potential 


and 


Q- 


in definitions of coordinates), 




r- 


df— 5.46, 




r e ~ 


D m — 5.44 (also 5.50 and 5.53), 




n,r 2 - 


e — 3.24; element of 2d set, 




Ri,R2,Rs 


e — 3.4; element of 1st set, 




Rm(z)- 


e'— 8.32, 




sgn ol t 


E — 5.12; "eccentric anomaly," 




t- 


E — 8.17; secular part of E( = M 8 ), 




v- 


E p — 8.17; periodic part of E, 




V 8 - 


Eq,E 1 ,E 2 — various terms of E p , 




Vv 


E(q) — complete elliptic integral of second kind, 


VQ,V h V 2 - 


modulus q, 




Va 


E(\p,q)- — incomplete elliptic integral of second 


kind; 


X,Y,Z 


6.9 and 6.10, 




x- 


/(p)— 2.10, 




y- 


F(p)— 2.8, 




a r 


F(\p,q) — incomplete elliptic integral of the first kind; 




6.9 and 6.10, 




a 2 - 


0(h)— 2.9, 




«3- 


h — 5.7, 






io — 3.5; element of 1st set, 




Pl,fo,fo 


/ — 4.7; element of 2d set, 




r 


I c — value of / for which X=l, 




7m" 


J 2 — 2.4; coefficient of second zonal harmonic of 


5 n " 


potential, 




€" 
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Symbol Definition or equation number 

rj — 2.1 and 2.2; a spheroidal coordinate — »sin 6 

as r— >co, 
7/0—3.33, 3.36, 3.39; during motion 
— 1 ^ — ?7o ^=77 ^770 = 1, 

772—3.33,3.37, 77 2 >>1, 
S — 2.1 and 2.2; geocentric declination, 
X— 5.6; bifa, 
*\,v2,vz — 7.4 through 7.6, 

p — 2.1 and 2.2; a spheroidal coordinate — >r 

as r— >oo ? 
<f> — 2.1 and 2.2; geocentric right ascension; the 
third spheroidal coordinate, 



b^mbol Definition or equation number 

x-6.51, 
\p — 6.4; a uniformizing variable analogous to 

the argument of latitude, 
\f/ B — 8.17; secular part of ^, 
\p v — 8.17; periodic part of \f/, 
^0,^1,^2— various terms of \p v , 

co — mean motion of p-perigee relative to the 
_ ascending node, 

i'l — mean motion of the ascending node relative 
to OX. 
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